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Abstract 

In 1952, Kac and Ward developed a combinatorial formulation for the two 
dimensional Ising model which is another method of obtaining Onsager's fa- 
mous formula for the free energy per site in the termodynamic limit of the 
model. Feynman gave an important contribution to this formulation conjec- 
turing a crucial mathematical relation which completed Kac and Ward ideas. 
In this paper, the method of Kac, Ward and Feynman for the free field Ising 
model in two dimensions is reviewed in a selfcontained way and Onsager's 
formula computed. 

Em 1952, Kac e Ward desenvolveram uma formulagao combinatorial do mod- 
elo de Ising em duas dimensoes que e um outro metodo para se obter a famosa 
formula de Onsager para a energia livre por sitio no limite termodinamico 
do modelo. Feynman fez importante contribuigao a esta formulagao conjec- 
turando uma relagao matematica crucial que completou as ideias de Kac e 
Ward. Neste trabalho, o metodo de Kac, Ward e Feynman para o modelo 
de Ising em duas dimensoes sem campo e revisada e a formula de Onsager e 
calculada. 
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1 Introduction. 



The aim of statistical physics is to understand the macroscopic behaviour of a system 
formed by a very large number of particles from information about how they interact 
with each other. One way in which one can gain insight into this problem and 
thus about complex systems is by constructing idealized models which hopefully 
will exhibit some of the interesting features of real systems like phase transitions. 
Perhaps the most studied of these idealized models is the Ising model so called in 
honor to his first investigator, Ernst Ising (1900-1998). 

The model was originally proposed as a simple model of ferromagnetism. In 
ref. [1] Ising investigated the model in one dimension and computed exactly its parti- 
tion function. In 1944, Onsager [2] considered the free field model in two dimensions 
and succeded to compute the partition function exactly. His method became known 
as the algebraic formulation of the model. In 1952, Kac and Ward [3] developed 
a quite different method of obtaining Onsager results known as the combinatorial 
formulation of the Ising model. Feynman developed the method farther and con- 
jectured an identity relating functions defined on graphs and functions defined on 
paths on a square lattice [4, 7]. This identity is a crucial element in the combi- 
natorial formulation of Kac, Ward and Feynman of the Ising model. The identity 
was later formally proved by Sherman [4-6] , followed later on by another proof by 
Burgoyne [7]. A somewhat similar treatment to the combinatorial formulation of 
Kac, Ward and Feynman can be found in refs. [12-14]. An important variant of the 
combinatorial formulation using the so called Pffafians was developed by Green and 
Hurst [10]. 

The bibhography on the Ising model is vast and to give a full list of references 
is virtually impossible. A nice introduction to the model though is the paper by B. 
Cipra given in ref. [17]. Old surveys but still useful on the distinct formulations 
of the Ising model in two dimensions and its history can be found in refs. [10-11, 
15-16] together with full lists of original references. 

The objective of the present paper is to review in a selfcontained way the 
calculation of the Onsager's formula for the two dimensional free field Ising model 
in the combinatorial formulation of Kac, Ward and Feynman. Our presentation 
follows chapter 5, section 5.4, of Feynman's book [9] and the paper by Burgoyne [7] 
although we have tried to be more careful with the mathematics involved than these 
references arc. 

The paper is organized as follows. In section 2, the Ising model is defined. 
In section 3 and through its various subsections the combinatorial formulation of 
Kac, Ward and Feynman of the partition function is given. In section 4, Onsager's 
formula for the free energy per site in the thermodynamic limit is computed. 
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2 Definition of the model. 



The model is defined on a finite planar square lattice A which mimic a regular 
arranjement of atoms in two dimensions. Suppose the lattice is embedded in the 
plane with sites having coordinates in Z x Z. To each site i of A it is assigned two 
possible states also called "spins" and denoted by (7,, where ai — +1 or ai — —1. 
The interaction energy between two particles located at the i-th and j-th sites and 
in the states CTj and aj, respectively, is postulated to be 




if i,j are n.n. 
otherwise 



(2.1) 

where "n.n" stands for nearest neighbors, hence, in the Ising model it is assumed 

that the energy depends only on short range interactions. The energy is —J if 
the nearest neighbors are in the same state and +J if the states are distinct. The 
constant J which can be positive or negative is a parameter for the model. 

Suppose A has A^^ sites. Then, there are 2^^ distinct configurations of the spins 
and, therefore, 2^ configurations a = (cxi, aN2) of the system. Call S = {a} the 
set of possible configurations of the system. The energy of each configuration a E S 
is given by 

E,^-J '^i'^j (2-2) 

n.n. Sit 

Suppose as well the system is at equilibrium temperature given by T. According to 
statistical mechanics, the probability to find the system in the configuration a is 

P, = ^f'^' (2.3) 
where /3 = ks is Boltzmann constant, and 

Z(/5)=^e-^^<^. (2.4) 

is the so called partition function of the model. This simple looking function is 
simple to compute exactly only in one dimension, difficult but possible to compute 
exactly in dimension two. In three dimensions nobody knows how to do it. 

The exact knowledge of Z{(3) allows one to obtain information about the global 
behaviour of the system. Important quantities that arc relevant to understand the 
physics of the system are all defined in terms of InZ or its derivatives. For instance, 
the free energy per lattice site / in the thermodynamic limit is defined as 

A basic problem is to find a closed form, analytic expression for /. Phase transitions 
will appear as singularities in / or in one of its derivatives. 
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3 The combinatorial formulation 



In the combinatorial formulation the partition function is expressed as a sum over 
special subsets of the lattice A called admissible graphs. Next, using a relation first 
conjectured by R. Feynman the resulting expression is converted into a product over 
paths. The final step towards the Onsager's formula to be accomplished in section 
4 consists in deriving an integral representation for this product. 



3.1 The partition function as a sum over graphs. 

Let's rewrite the partition function (2.4) as 

zn{k) = e • • • e n ^"""'^^ (3-1) 

CTl=±l CTJV=±1'»-"- 

with K — +]^- Noting that OiUj — rkl, it follows that 
and 

JJe^'^*'^^- = (l-ii^)-i \{{l + a,<jju). (3.3) 

n.n. n.n. 

where u = tanhi^ and x = 2N{N — 1) is the number of bonds in A. Notice that 
I M |< 1, for any K. 

Definition 3.1. An admissible graph is a connected or disconnected subset of A 
whose sites have even valence. 

Definition 3.2. Given an admissible graph G, define 

Ig{u)^IIu^u'^ (3.4) 

where the product is over the bonds i of G. 

Theorem 3.1. Call A the set of all admissible graphs G of A. Then, 

Z^{u) = 2^\l - uY^^^-'^ (l + E W] (3.5) 

V GeA I 

Proof: To each pair i,j of nearest neighbors of A there correspond a term uaiaj 
and a bond. Since the number of pairs i,j of n.n. coincide with the number x ~ 
2N{N — 1) of bonds of A the product on the RHS of (3.3) is a polinomial of degree 
X, that is, 

X 

l[{l + uaiaj) = l + E"^E(^n^i2)---((^i2p-iC^i2p) (3-6) 

n.n. p=l n.n. 
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The second summation is over all possible products of p pairs {(TiCTj) of n.n. of 
A where a pair is not to occur twice in the same product. To each pair 
there is associated a bond connecting the neighbors i and j so to each product of p 
pairs correspond a graph (connected or disconnected). So, the second summation is 
over all graphs with p bonds. The graphs may have sites with valence 1,2,3 or 4. 
The summations over the spins (Tj's eliminate graphs having sites with odd valence 
because "^j = ^-nd J2 o'f = 0. The graphs left arc those whose sites have valence 
2 or 4, thus admissible. If Vq is the number of sites in a admissible graph G then 
there is a factor 2^° associated to it because each site of G contributes a factor 2 
coming from — 2 and J^crf — 2. The sum over a includes all the (Tj and not 
only those with sites i in G. The summation over the sites not in G will give a factor 
2^~^G where V — N"^ is the number of sites in the lattice, hence, in the end one gets 
the factor 2^. 



3.2 The partition function as a product over paths. 

Let's orient and number the bonds of A with distinct positive integers i and call A 
with this indexation a colored lattice. 

Definition 3.3. A path p over A is an ordered sequence of bonds each starting at 
the site where the previous one ended. The last bond ending at the site from which 
the first one started. Thus, p is closed. The path is subjected to the constraint that 
it never goes backwards over the previous bond. A path p is given by a word, that 
is, an ordered sequence of symbols Dj where i distiguishes the bonds of A. A path 
p is then of the form 

P-D^lD^l..D]l (3.7) 

for some I and where Cj = +1 (— 1) if the path traverses bond ji following the direction 
(opposite direction) assigned to it. Because a path is closed it is defined to within 
its circular order so that 

DjlD^l...D]l ^ D^...DID'^ ^ ... = DlD-l...iy^ (3.8) 

The inversion p~^ of p is given by 

DTD -Dr (3-9) 

We take p and p~^ to be equivalent. Given p, denote by [p] the set of all paths 
equivalent to p, that is, its circular permutations and their inversions. 

Definition 3.4. A periodic path is one which has the word representation 

{D]l...D]ir (3.10) 
for some I and w >2 and where the subword in between brackets is nonperiodic. 
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Definition 3.5. A path p has assigned to it a sign given by 

s{p) = (-1)^+* (3.11) 

where t is the number of 27r-angles turned by a tangent vector while traversing p. 
A positive (negative) angle is assigned to a counterclockwise (clockwise) rotation. 

Example 1. See Figure la). A tangent vector starting at point e and traversing 
the path shown in Figure la) turns once a total angle given by 4.| = 2tt after its 
return to e so in this case t = 1 and s{p) = +1. For the path in Figure lb), the 
total angle turned is 3| — 3| = so i = and s{p) — —1. 



e e 
(a) (b) 

Figure 1: Examples of paths with (a) s{p) = +1 and (b) s{p) = —1. 

Renicirk. In section 4, instead of assigning an angle ±7r/2 to a turn we will count 
the contribution to the sign by assigning a = e"^^^^ and a = e"""^/^ to each coun- 
terclockwise and clockwise turn, respectively, and then in the end multiplying the 
result by —1. In the example above, one gets in this manner (e^'"^^)^ = —1 and 
^gi7r/4^3^g-j7r/4^3 _ _|_-|^ Multiplying both results by —1, one recovers the correct 
sign for each path. 

The sign of periodic paths. Suppose the sign of the nonperiodic path in between 
brackets in (3.10) is (— 1)-'^+*^ Then, the sign of the periodic path with period w is 

(— 1)^+"'*^ Hence, the sign of a periodic path is —1 if its period is an even number 
and the sign equals the sign of the nonperiodic subpath if the period w is an odd 
number. 

Definition 3.6. To each path p it is assigned the function Ip{u) given by 

Ip{u) = (3.12) 

where / = mi + ... + ruk, for some k, is the length of p, rrii being the number of times 
bond i is covered by p, and the function Wp{u), "the amplitude of p", defined as 
follows: 

Wp{u) = s{p)Ip{u) (3.13) 
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Theorem 3.2. The functions laiu) and Wp{u), \u\ < 1, defined above satisfy the 
foUowing relation: 

1+ ^/GH=n[i+^pH] (3-14) 

GeA [p] 

The product is over all inequivalent classes [p] of closed nonperiodic paths. The 
summation is over all admissible graphs of the finite N x N planar square lattice A. 

Relation (3.14) is a simpler version suitable for the Ising model of a more 
general relation investigated by Sherman and Burgoyne in refs. [4-7] . The difference 
is that they assign to the bonds i of the lattice distinct parameters di, hence, in this 
case the functions Iq and W are given in terms of these parameters. In the Ising 
model context under consideration these are all equal to u and \u\ < 1. 

According to references [4,7,10,11], relation (3.14) first appeared as a conjec- 
ture in lecture notes by Feynman ( ref. [9], published only in 1972 and already 
mentioning ref. [4]). The first proof of it was achieved by Sherman in refs. [4,6] 
followed by another one later on by Burgoyne in ref. [7]. The simplest nontrivial 
case of the general relation is investigated in ref. [8]. 

Below Burgoyne's proof is essencially reproduced for the case \u\ < 1. 

Proof: Expand the product over the distinct classes of nonperiodic paths [p] as 1 
(one) plus an infinite sum of terms of the form 

Wp,{u)Wp,{u)...Wp,{u) ^sl[u'' (3.15) 

i 

for some k where pi, is a set of nonperiodic paths over A. The product on the 
r.h.s of (3.15) is over the bonds i traversed by pi,p2, ...pk, and says how many 
times. Ifpi,p2, ■■■Pk traverse bond i, say, mi{i), ■■■,mk{i) times, rrij > 0, respectively, 
then = Yli=i iTT'ji'^)- The sign s is the product of the signs of pi,P2, ■■■,Pk- 

Let's prove, first, that those terms having rj = 1, Vi, add up to Z]/g(m). 
Consider one of these terms with associated paths Pi,P2, ■■■,Pk- Each bond in the 
set of bonds traversed by Pi.p2, ■■■iPk is traversed only once by one of these paths. 
Thus, the only possible intersection if any between any two of these paths in this 
case can occur only at a site of valence 4 and they cross each other like in Fig. (2. a) 
. Otherwise, they are disjoint. Thus, the set of bonds traversed by paths Pi, 
constitute a graph whose vertices have valence 2 or 4. This is an admissible graph. 
Therefore, to each term of the form of (3.15) with = 1, Vi, one can associate 
an admissible graph. This graph can be disconnected. This happen if the set of 
paths can be split into subsets completely disjoint which generate admissible graphs 
without any bonds and vertices in common. 

Now, given an admissible graph G one can in general associate more than one 
term of the form of (3.15) with rj = 1, each associated with a distinct set of paths. 
Let's see how this follows. The sites of an admissible graph have valence 2 or 4. 
When a path strikes a site of valence 4 it has only 3 possible directions to follow. 
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See Figures 2a, 2b, 2c. (The case in Fig. 2d is forbidden.) Then, any two terms 
associated to a given admissible graph G wiU differ only in the types of crossings 
at the sites of G. Since there are 3 types of crossing per valence 4 site, the number 
of possible terms associated to G is 3^ where V is the number of sites of G with 
valence 4. 

A term has a sign which comes out from the contribution of the signs of the 
paths associated to that term. Let's see how the sign of a term comes out. A term 
with ti crossings of type 1 (Fig. 2a.) has a sign which can be expressed as (— 1)*^ 
where ti includes selfcrossings of single paths plus crossings between different paths. 
Indeed, since distinct closed paths always intersect in a even number of crossings 
then (—1)*^ will give the correct sign of the term which is the product of the signs 
of the individual paths. Let's associate to the crossings of type j = 2, 3 the sign 
(+1)*2(+1)*3 so that a term with tj crossings of type j — 1, 2, 3 has a sign given by 



There are V\ ways of distributing V = ti + 12 + ^3 crossings among the sites of 
G but since there are tj crossings of the type j, j = 1, 2, 3, one has to divide V\ by 
ti\t2lt3\ SO that the number of distinct terms with tj crossings of type j is 



These terms have the same factor I{G) — where L is the number of bonds of G. 
Summing all these terms arising from a given G and summing over all admissible 
graphs G of A the result is 



where J^t means summation over all ti, ^2, ^3 such that ^1 + ^2 + ^3 = V. Using the 
multinomial theorem the summation over {t} gives (—1 + 1 + 1)^ and one gets the 
result Eg^(G'). 

If G is disconnected with I components Gj, i — 1, 2, Z, each of them with tj, 
j — 1, 2, I, sites of valence 4 and J^tj = V, then applying the previous argument 
to each component will give I {Gi)I {G2) {Gi) — I{G). 

In view of the above result, the theorem could be equivalently stated by saying 
that the sum of terms with Tj > 1 for at least one of the i converges to zero. Let's 
prove this. 

Let Q be the set of all colored connected or disconnected subgraphs g of the 
colored lattice without valence 1 sites and such that if g is connected then g is not a 
poligon, that is, a graph having valence 2 sites only. A disconnected graph is allowed 
to have some but not all of the components as poligons. The reason for excluding 
graphs which are poligons or having all components which are is that closed paths 



(_l)*l( + l)t2(+l)t.. 




(3.17) 
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with repeated bonds over them are necessarily periodic and these are forbidden. 
The coloring of g is that inherited from the colored lattice. 

Given g & Q, call ii, ii(^g) the bonds of g. A term Wg associated to g is of the 

form 

^9 = ^Pi^P2--^Pk = {signwg)\wg\ (3.18) 
for some k and set of paths pi, ...,pk, which traverse the bonds of g only, where 

Kg) 

\wg\ = l[u'^ (3.19) 

and is the number of times bond ij is traversed by pi, that is, li pi, ...,Pk 

traverse the i-th bond mi{i), ...,mk{i) times, m > 0, respectively, then 

k 
a=l 

Some but not all of the m's can be zero so that rij > 1 with at least one > 1. 

Let's consider the set of all terms with the same effective set of bonds {i} 
traversed, hence, the terms associated to a given g. Within this set it's possible in 
general to find terms with the same powers {r} and the same \wg\ although having 
distinct associated paths and possibly with different effective sign. 

Let's group together those terms which cover the same bonds of g the same 
number of times. Denote by VVg,Ar(r) the set of terms Wg with the same powers 
{vij } and such that X]j=i ^ij = for fixed N. The summation over all terms with 
repeated lines can now be expressed as 

E E E E (3-21) 

geg {N}g r{N) WgeWg,N{r) 

where J2geg i^icans summation over all elements in Q; Y^{N}g means summation over 
all positive integers compatible to the given graph g and such that N > l{g) + 1; 
J2r(N) means summation over a set of positive integers ri, such that ri + ...+r/ = 
N and which are also compatible to g; and, finally, J2wg means summation over all 
terms Wg e Wg,iv(0- 

Now the following remarks come to order. In the second summation, the case 
N — I is excluded for it implies that ri — 1 and in this case there can be no repeated 
bonds. The case N < I corresponds to another element g' e Q. The equality depends 
on the graph g. For instance, take the graph shown in Fig. lb where / + 1 =9. 
No nonperiodic closed path with repeated bonds can have length = 9 because 
l{g) = 8. The length N can only be even and its minimum is N — 12. Hence, for 
this particular graph g the summation is over all even numbers greater or equal to 
12. In any case, the set {N}g has always infinite elements. Given g and N e {N}g, 
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not all partitions of N are allowed in the third sum. For instance, given the graph 
in Fig. lb and N = 12, the partition with r^^, = 1, VA; ^ 1, and rj = 5 can not 
be associated to any allowed path. So, the set of integers {N} and partitions of N 
must be suitable to each g. 

Given g, let's consider now the partial sums 

E E E (3-22) 

{N\N<n}g r{N) WgeWg,N(r) 

The goal is to show that in the limit n — > oo, s„ goes to zero. In ref. [7] its proved 
that Sn — 0. The argument of the proof goes as follows. 

Since the bonds of g are covered the same number of times by all elements in 
the group, choose a bond of g, say b, which is traversed f > 1 times by all elements 
in the group Wg^N- This choice has to be done for each partition r{N). Denote 
by P the set of paths associated to Wg. Then, P = P'[jP" where P' is the set of 
those paths which traverse bond b whereas P" is the set of those paths which do not 
traverse b. 

Given a path p E P', let Pc be the path segment obtained from p upon removal 
of b. Given P' — {p,p'...} define 

Pc = {Pc,p'c,-} (3.23) 
This set has exactly f path segments. Collect under a same subgroup S the elements 
Wg e Wg,N having the property that line b is covered exactly f times by all elements 
in S and they all have the same subset P^ with r path segments and the same subset 
P". The set Wg^N is the union of such subsets, that is, 

E E E^HKI (3-24) 

where s{w) is the sign of Wg and \wg\ = . Recall that |m| < 1 so that |-u;g| < 1. 

The elements inside any given S cancel each other. Denote by q and e the 
elements of S that are in P' and P", respectively. Suppose that the segments 
Qi, ...,qr are all distinct. (For the case with repeated segments, see [7]). The terms 
in S are precisely those which can be obtained by joining the ends of the segments 
and this can be done in exactly r! ways. This gives the possible terms Wg in the 
subgroup. From the properties of the permutation group half of A^! permutations 
are odd and half are even and so the signs of half of the terms are positive and half 
are negative, hence, a cancellation takes place. 

Using (3.14), the partition function of the two dimensional Ising model can 
now be expressed as a product over paths as follows: 

Z^{u) = 2^'(1 - «2)-iv(iv-i) -Q ^ ^3 25) 

[p] 

The next step consists in expressing the product over [p] as an integral. This will 
be achieved in the next section. 
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4 Paths amplitudes and Onsager's formula 



Consider all paths that start at a fixed site Pi which we take as the origin with 
coordinates (0, 0) and end at the site Pn+i with coordinates {x, y) in n steps. Starting 
at (0, 0) and whenever a site is reached there are four possible directions which a 
path can take (see Figure 2 and the Remark below) . The path a) continues forward 
in the same direction of the previous step; b) it turns left 90" relative to the previous 
step; c) it turns right 90'^ relative to the previous step; d) it turns 180°. To each 
one of this possibilities it is assigned an "amplitude" which is given by: A) u for 
the case a); B) ua for the case b); C) ua for the case c) and D) (zero) for the 
case d) , where u — tghK and a — e*'^/^ is the contribution to the sign of p each 
time it turns left (counterclockwisely) relative to the previous step and a, when it 
turns right (clockwisely) . See the Remark after Example 1, sec. 3.2. 



(a) (b) (c) (d) 

Figure 2: Directions which a path can take at a valence 4 site. 

Remark. The lattice being finite it has a border so that when a path strikes a site 
on the border it may have there only two or three possible directions to follow. In 
the spirit of refs. [7,9] we shall neglect the border and derive the relevant formulas as 
if there was no border at all with the justification that in the limit — > oo which we 
shall take in the end of the calculation border effects dissapear. Of course, another 
approach would be to do everything on a toroidal lattice. In this case, however, 
relation (3.14) must be replaced by another more involved identity apropriate for 
the toroidal lattice ( given in refs. [4, 10] ). We shall restrict the presentation to 
the planar case only. 

Call Un{x,y) the amplitude of arrival at {x,y) moving upward in the n-th 
step, Dn{x, y) the amplitude of arrival at (x, y) moving downwards in the n-th step, 
Ln{x, y) the amplitude of arrival at (x, y) moving from the left in the n-th step, and 
Rn{x, y) the amplitude of arrival at (x, y) moving from right in the n-th step. 

If the path arrives at (x, y) moving upward in the n-th step then 

Unix,y) = uUn-iix,y - 1) + ODr,-i{x,y - 1) +uaLr,-i{x,y - 1) +uaRn-i{x,y - 1) 

(4.1) 

where Un-i, D^-i, L^-i and Rn-i are the amplitudes associated to the four pos- 
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sibilities to reach site {x,y — 1). Relation (4.1) can be understood as follows. If 
{x,y — 1) is reached going up a bond in the (n — l)-th step, there the amplitude 
is Un-i{x,y — 1) so in the n-th step as the path follows the same direction of the 
previous step, by the rules a) and A) above, a factor u is multiplied to the amplitude 
Un-i{x,y — 1). See Figure 3. 

T (x,y) 



(x,y-i) 



Figure 3: p goes up to (x, y — 1) and (a;, y) in the (n — l)-th and n-th steps. 



If the site (x, y — l) is reached from the left in the (n — l)-th step (Figure 4), the path 
has to make a counterclockwise rotation to go to (x, y) in the n-th step. By the rules 
h) and B) a factor ua should then be multiplied to the amplitude L„_i(x, y — l). 



> i (x,y-l) 

Figure 4. p turns counterclockwise at (x, y — 1) to go to (x, y) in the n-th step. 



The case that the path goes down to (x, y — 1) in the (n — l)-th step and goes up 

to {x,y) in the n-th step corresponds to a 180*^ rotation. By rules d) and D) the 
amplitude should be 0.-D„_i(,t, y — l). If the site (,t, y — 1) is reached from the right 
the path has to make a clockwise rotation to go to (a;,?/) (Figure 5). By the rules 
c) and C) a factor ua. should then be multiplied to the amplitude Rn-i{x, y — l). 
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(x,y) 



(x,y-i) 

Figure 5. p turns clockwise at {x, y — 1) to go to {x, y) in the n-th step. 

Analogously, if a path arrives at {x, y) in the n-th step going down the amplitude is 
given by the relation 

Dn{x, y) = OC/n_i (x, y + 1) + uDn-i (x, y + 1) + uaLn-i (x, y + 1) + uaRn-i{x, y + 1) 

(4.2) 

If it arrives at (x, y) coming from the left then the amplitude is given by 
Ln{x, y) = uaUn-i{x - 1, y) + uaDn-i{x - 1, y) + uL^-iix - 1, y) + ^Rn-i{x - 1, y) 

(4.3) 

At last, if it arrives at (x, y) coming from the right the amplitude is 

Rn{x,y) = uaUn-i{x + l,y) +uaDn-i{x + l,y) + OLn-i{x + l,y) +uRn-i{x + l,y) 

(4.4) 

Of course to compute an amplitude using the above recursion relations it is needed 
the amplitude in the zero-th step. We shall follow the convention of reference [10], 
namely, that in the zero-th step a path arrives at the origin moving upward so that 
Uo{x,y) = SxflSy^o and Dq = Rq = Lq = 0. The amplitude to arrive in zero steps is 
one if the path arrives going upward at the origin and zero for any other point or 
any other direction of arrival. 

Example 2. See Figure 6. Let's compute the amphtude of arrival at site (2,1) in 3 
steps moving upward in the third step. Only one path is possible in this case. Using 
the recursion (4.1), 

C/3(2, 1) = uU2{2, 0) + 0D2(2, 0) + uaL2(2, 0) + uaR2(2, 0) (4.5) 
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(2,1) 



(0,0) (1,0) (2,0) 



Figure 6: The path in Ex. 2. 



In the second step, the path moves to site (2,0) coming from the left so U2 — D2 — 
i?2 = and C/3(2, 1) = uaL2{2, 0). From (4.3), 



with [/i Di = i?i = so that U^{2, 1) = u^aLi{l, 0) where 

Li(l, 0) = uaUo{Q, 0) + iiQ;Do(0, 0) + wLo(0, 0) + 0i?o(0, 0) = ua (4.7) 
implying that U^{2, 1) = v?. 

Example 3. Let's now compute the amplitude of arrival at (2,1) in 3 steps moving 
from the left in the third step. In this case, the possible paths are shown in Figure 
7a) and 7b). 



L2(2, 0) = u6iUi{l, 0) + uaDiil, 0) + uLi{l, 0) + Oi?i(l, 0) 



(4.6) 



(1,1) (2,1) 



(0,1) (1,1) (2,1) 



(0,0) (1,0) 



(0,0) 



(a) 



Figure 7: The paths in Ex. 3. 



Using relation (4.3), the amplitude is 



L3(2, 1) = uaU2{l, 1) + uaD2{l, 1) + 1^^2(1, 1) + Oi?2(l, 1) 



(4.8) 



Using (4.1), 



C/2(l, 1) = uUi{l, 0) + ODi(l, 0) + uaLi{l, 0) + uaRi{l, 0) 



(4.9) 



14 



Since Ui = Di = Ri = 0, one finds that C/2(l, 1) = uaLi{l, 0) = uaua — v?. Using 
(4.3), with Di = Li = i?i = 0, 

L2(l, 1) = u6iUi{Q, 1) + iiQ;Di(0, 1) + 'uLi(0, 1) + Oi?i(0, 1) = uau (4.10) 

Therefore, L^{2, 1) = 2u^a. 

Definition 4.1. The partial amphtude of a path p of length n is given by 

Wp{u)^J{aJ{au'' (4.11) 

Definition 4.2 . The amplitude Z^p W^p(n,Pi)('w) of arrival at P„+i(a;, y) from any 
direction in n steps is given by 

Un{x,y) + Dn{x,y) + Ln{x,y) + Rn{x,y) (4.12) 

Example 4. The partial amplitudes for the paths in Figure 6, 7a) and 7b) are m^, 
au^ and aii^, respectively. The amplitude of arrival at (2,1) from any direction in 3 
steps is, then, + 26iu^. 

Definition 4.3. Fix n and call C„(x, y) the set of all paths of length n starting at 
(0, 0) and arriving at (x, y). Given p & Cn and F„ e B{x, y) where 

^{x.y) = {Unix,y),Dnix,y),Lnix,y),Rnix,y)} (4.13) 

Define the extension of Fn{x. y), denoted by the same symbol, so as to include sites 
(x, y) which can be reached only by a number m > n of steps but in this case set 
F„(x,y)=0. 

Lemma. The transform of F„, the function F„(e, 77), < e < 27r and < 77 < 27r, 
given by 

00 00 

^n(^,r7)= E Pn{x,y)e-^^^e-'^^ (4.14) 

x=—oo y=—oo 

is well defined and 

F„(x, y) = fj fj e-^e^''^F„(£, ry)^ (4.15) 

Proof: Fn{x,y) = for | x |> n or/and \ y \> n. Then, for fixed n the sums in 
(4.14) have only a finite number of terms. 

Using (4.14), the transform of Un{x,y) is: 

00 00 

Un{s,rj)^ Y: E t/n(a;,y)e--e-''^ (4.16) 

1=— 00 y=—oo 

Upon substitution of (4.1), and making the change y = y — 1 ii follows that 
[/„(£, 77) = we""'C/n_i(£, 77) + OL>n_i(£, 77) + uae'^"^ Ln-i{e , 77) + uae~^'^Rn-i{e, 77) 
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(4.17) 

Similarly, we obtain Dn{£,r)), Ln{s,r)) and Rn{s,r)): 

(4.18) 

Ln{e, rj) = uae~^^Un-i{e, r]) + uae'^-^D n-i{e, rj) + we~*^L„_i(£, rf) + QRn-i{e, 77) 

(4.19) 

'Rn{e, 7]) = uae'^Un-i{e, 77) + uae'Wn_i{e, t]) + OL„_i(£, r/) + ue'^Rn-iie, rj) 

(4.20) 

Call ipn{^,v) the matrix 

Un Dn Lfi Rn ) (4.21) 
Then, from (4.17-20) we obtain that 

V'n(£, V) = ^n-i{e, r])uM, (4.22) 

where 



M = 



f V ah ah \ 

V ah ah 

av av h 

\ av av /i y 



(4.23) 



with V = e h = e h = e**^, v = e*'' and a = e** . 

Call 1,2,3 and 4 the directions shown in the Figure 8 below: 



Figura 8. Directions associated to Mj 



Notice that the subindices i,j of M^j are in one-to-one with the directions. Indeed, 
uMij corresponds to the amplitude of arrival at {x, y) ( in (e, rf) space) coming up 
in the {n — l)-th step, Vj, but going up if j = 1, down if j = 2, coming from the 
left if J = 3 and coming from the right if j = 4 in the n-th step. Therefore, uMij 
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is the amplitude of arrival at (x, y) ( in (e, rf) space ) following directions 1 and j in 
the (n — l)-th and n-th steps, respectively. More generally, uMij is the amplitude 
of arrival at (a;, y) following directions i and j in the (n — l)-th and n-th steps, 
respectively. 



Prom now on only closed paths starting at (0,0) and arriving at (0,0) in n steps 
will be considered. 

Prom (4.22) it follows that 

V'n(£,77) = 'il^n-i{e,r)){uM) = 'il;n-2{e,r)){uMf = • • • = Vo(^iM)" (4.24) 

Denote by ■0o,j, 1 < ^ < 4, the line matrix with the only element distinct from zero 
and equal to 1 in the i-th column. Let ipQ = ijjo^i according to whether the path 
arrives at the origin moving up (i = 1), down (i = 2), from the left {i = 3) or from 
the right (i = 4), respectively. Then 

F^{e,ri)=^o,i{^Mri^l, (4.25) 

where i = 1, 2, 3, 4 if F = [/, D, L, R, respectively, and is the transpose of ^ and 

Fn{e,v) = Ei'oA^,v){uMri'o,kie,v) (4-26) 
F„eB„ fe=i 

Given a 4 x 4 matrix A, ipo^iA is the line matrix formed by the elements in the i-th 
hne of A, that is, 

jjjQ^iA = ( Ai^i Ai^2 Ai,3 Ai^i ) (4.27) 
so ipo^iAip'^^- — An. Therefore, the sum over i equals the trace of A. Thus, 

^ ^lJo,k{uMY'il^lk = Tr{uMr (4.28) 

k=l 

The total partial amplitude of arrival at (0, 0) of closed paths moving in any direction 
in n steps given by (4.12) can be expressed compactly as 

E (4.29) 

Fn&Bn 

Prom (4.15), (4.26) and (4.28), it follows that 



To better understand relation (4.30), consider the matrix u"'M^, for some n. 
An element {u^M"')ij^i^_^_-^ of this matrix is given as 

4 

(u^M^),^,^^, = E uM,^,,uM,,,, . . . mM,„_,,„mM,„,„^,. (4.31) 
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Recall that uMij is the partial amplitude of a path arriving at a site coming from 
direction i and going to the next site in one step following direction j. Thus, each 
term in the r.h.s. of (4.31) is the amplitude of a path of length n starting at Pi 
coming from direction ii, going to P2 following direction 12, etc, and arriving at site 
Pn+i following direction in+i after n steps. The element {u^M'^)i^i^_^^ gives the total 
partial amplitude of arrival at Pn+i in n steps in (e, 77) space. 

The terms in (ii"M")jjj^_^j describe open as well as closed paths. Let's see 
some examples. 

Example 5. Take n = 5, ii = 2 and — The term M23M31M11M14M41 describes 
a path beginning at Pi where it arrived coming from direction ii = 2, going to P2, 
P3, P4, P5 and to Pq following directions ^2 = 3, is = 1, ^4 = 1, ^5 = 4 and = 1, 
respectively. See Figure 9a) . 



Pfi 



Pi 



Pa 



P. 



Pi 



Pa 



(a) (b) 
Figure 9: Paths in (a) Ex. 5 and in (b) Ex. 6. 

Example 6. Take n = 6, ii = ie = 2 and the term M23M31M11M14M42M22 of 
(M^)22 . This term describes the closed path in Fig. 9b. The elements of 
outside the diagonal have associated to them only open paths. This is implied by 
the simple fact that these elements have ii ^ in+i- Closed paths are to be found 
only in the diagonal elements since there ii = in+i- However, open paths can also 
be associated to some terms in the diagonal elements. Let's see some examples. 



Example 7. Take n = 2, ii 
M13M31 + M14M41 with 



i3 = 1 and the element (M^)ii = M11M11 + M12M21 + 



'u^(M^)ii = + + u^{ah)(av) + u'^(ah)(av) 



(4.32) 



To each one of the terms of (M^)ii correspond the paths (a), (b), (c) and (d), 
respectively, shown in Figure 10. 
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P3 
Pi 



Pi 



-Pi Pi - 
t t 



P. 



(a) (b) (c) 

Figure 10: Paths associated to the terms in (M^)ii 



(d) 



Example 8. Take n — A and consider the following terms in (M^)ii: 

a)The term u^MnMi-iMuMu — u^v'^ is the amplitude of the open path shown in 
fig. 11 below. 



Figure 11: Path (Mn)! 

b) The term M11M11M13M31 = vv{ah){av) — v^hv whose associated open path is 
shown in Figure 12. 



Ps 

P2 



Pi 

Figure 12: Path MnMnMigMgi. 

c) The term u^Mi3Ms2M24M4i = u^{ah){av){ah){av) = u'^a\hh){vv) 
the amplitude of the closed path shown in Fig. 13: 



4 — 4 • 
u a IS 
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Figure 13: Path M13M32M24M41. 

In order to restrict to the elements of Af " having closed paths we must take 
the trace of M". A closed path begins at and return to Pi after n steps. Since it is 
closed it has to cover n/2 horizontal bonds in one direction and n/2 horizontal bonds 
following the opposite direction. The same is true for the vertical bonds traversed 
by p. So, if the term Mi^i^Mi^i^ . . . Mi^_^i^Mi^i^^^ , in+i = in, describes a closed 
path, then the number of /I's 's) equals the number of /i's (u's) appearing in it. 
In this case, it's possible to organize the term into a product of pairs hh = 1 and 
vv = 1 and the double integral in e and 77 will give (27r)^ times a product of a's and 
a's. More precisely, the double integral over a closed path in TrM" equals 

(271)^1] "H^ (4-33) 

where the first product is over all counterclockwise rotations and the second is over 
all clockwise rotations, so 

l[aY[a = {-ly^P^ (4.34) 

where t{p) is the number of complete 27r revolutions performed by a tangent vector 
traversing the closed path p. Remind that one has yet to multiply (4.34) by (— 1) 
in order to get the complete sign s{p) of p. 

If a path is open the /i's, /i's (v's, v's) don't match up into pairs. There will 
be left integrals of the form 

f " e'^'^de = 0, (4.35) 
Jo 

where 6 stands for 77 or e and k > 1, hence, the integrals in rj and e remove completely 
terms describing open paths. Let's see examples. 

Example 9. Take n — 1. In this case there are only open paths and 

J J dr]d£Mij = 0, (4.36) 
Example 10. Using (4.35), in ex. 7, 

/ / dr]de{M^)n = (4.37) 

Jo JQ 

Example 11. It is clear that 

J dedr]{M'')ij ^ (4.38) 
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if n = 1, 2, 3, which is guaranteed by the fact that in a square lattice closed paths 
are possible only if n > 4. In the case n = 2 the path in Figure 10. b) is closed but 
it traverses the same edge back and its amplitude is thus zero. 

Example 12. Using (4.35) in ex.8, for the term MuMuMuMu = 



r27V r27T pZTT 

/ / v^dride = 2tt / e'^'^dr] - 0. 
Jo Jo Jo 

For the term M11M11M13M31 = v% 

/•2vr 

/ drjdev h = 
Jo 

For the term M13M32M24M41 = which describes a closed path the result 

1 r2iT 



drjdea 



(4.39) 



(4.40) 



(4.41) 



(27r)2 Jo 

follows which has the form (4.33-34) with = 1. 

Given a closed path in {M^)u, the inverse path is present in some {M")jj, 
j 7^ i. For instance, in (M^)ii there are the closed paths shown in Fig. 14 given by 
the terms M14M42M23M31 and M13M32M24M41. 



Pi 



Pi 



(b) 



Figure 14: Paths (a) M14M42M23M31 and (b) M13M32M24M41. 



In (M^)22 there are the terms M24M41M13M32 and M23M31M14M42 , with 
associated closed paths shown in Fig. 15c and 15d, respectively: 



Pi 



Pi 



(c) 



(d) 



Figura 15: Paths (c) M24M41M13M32 and (d) M23M31M14M42 
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In (Af^)33, there are the terms AL^2M24M4iMri and Af3iMi4M42M23 with as- 
sociated closed paths shown in Figure 16e) and 16f), respectively. 



'Pi 



Pi 



(e) 



Figura 16. Paths (e) M32M24M41M13 and (f) M31M14M42M23 



In (M^)44, there are the terms M42M23M31M14 and M41M13M32M24 with the 
associated closed paths shown in Figure 17g) and 17h), respectively. 



Pi 



Pi 



(g) 



(h) 



Figura 17. Paths (g) M42M23M31M14 and (h) M41M13M32M24 



Note that (e) is the inversion of (a), (f) is the inversion of (c), (g) of (b), and 
(h) of (d). 

So restricting to the diagonal terms of M" which amounts to take the trace 
of this matrix and then performing a double integration on the angles to eliminate 
open paths, dividing the result by 2 to eliminate inversions, and multiplying the 
result by —u"' gives the total complete amplitude ( with the right signs ) to arrive 
back at Pi in n steps moving in any direction. We have thus achieved the following 
relation: 

Y: W,{u) = -^7^ r r dedrjTr{uMr (4.42) 

p(„,PO 2(27r)2 7o Jo 

The above result is restricted to a fixed site Pi. For the finite N x N lattice with A^^ 
sites and disregarding boundary effects, the total (independent of site) amplitude of 
closed paths of length n is: 

1 /•27r /•27r 

J: W,{u)^J:WA^)--^77^2 / dedr]Tr{uMr (4.43) 
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Taking all A^^ sites into account imply that given a closed path p{n), the summation 
J2p(n) W{p) includes all circular permutations of p. To eliminate these the previous 
relation has to be divided by n. Then, the amplitude is given by 

Wj,{u) = - — 77^ / / dedr] \ (4.44) 



2 (27r)2 Jo Jo n 

We notice that a nonperiodic path appears n times in the sum but a periodic path 
of length n and period w has n/w distinct starting points only and for this rea- 
son it appears n/w times in the sum over paths. For instance, the periodic path 
{D j-^D j^){D j^D j,^){D j-^D j^) of length n = Q and period -u; = 3 has only two distinct 
starting points. The other equivalent periodic path is Dj^{Dj^Dj^)[Dj^Dj^)Dj-^^. Af- 
ter division by n, periodic paths with period w will show up in the sum with a weight 
1/w. Thus, the above relation includes all closed paths of length n over the N x N 
lattice, periodic and nonperiodic, and excludes inversions and circular permutations. 
The total amplitude of closed paths of any length is then given by the series 

E - E = E -^77^ / / dedri \ (4.45) 



^1 2 (27r)2Jo io n 

whose convergence will be investigated below. We note that since the lattice is 
square, closed paths with nonzero amplitude are possible only for n > 3 but in view 
of relation (4.38) in Ex. 11 we can write the series in (4.45) starting from n — 1. 

With the above remarks, 

00-1 -1-1 

E - E ^ UWpiu) - + i^iwMf - -] (4-46) 

n=l Pin) [P] 

In J2\p] the first term is the sum of Wp{u) over all nonperiodic paths. The other 
terms give the sum over all periodic paths since any periodic path is the repetition 
of some nonperiodic path p with period given hj w = 2,3,.... In section 3.2 the 
sign of a periodic path was proved to be —1 if w is even and equal to the sign of its 
nonperiodic subpath if w is odd. This explains the signs in the r.h.s of (4.46). 

Since \u\ < 1 then |II/| < 1 and the series between brackets converges to 
ln{l + W), and the r.h.s. of (4.46) equals to 

^ln(l + Wp{u)) = lnn(l + Wp{u)). (4.47) 
p p 

a result to be used below. 

Theorem 4.1. Take \u\ < r < j. Then, the series 



f (4.48) 

n=l 
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converges uniformily. 

Proof: We have that \Mij\<l, Vi, j = 1, 2, 3, 4, so from (4.31) we get 

I (M")ij |<4"-^ (4.49) 



and 



n 



4"-l I 7/ 1" 

< ^ (4.50) 



n 



The series 



n=l ^ 



converges for \u\ < r < j, hence, by Weierstrass M-test the series (4.48) converges 
uniformly for |m| < r < 1/4. 

We may conclude that the series Z^^i converges uniformly to the ma- 

trix ln{l — uM) in the same interval. 

We may now integrate the series term by term to get the series (4.45) which 
likewise converges uniformly in the same interval. Interchanging integration and 
summation in (4.45) yields 



-ITTTTT^ / dedrjTry ^ ——— / dedriTrlnil - uM) 

2 (27r)2 Jo Jo n 2 (27r)2 Jo Jo ' ^ ' 

(4.52) 

Prom the previous analysis, |m| < r < 1/4. However, the r.h.s. of (4.52) is well 
defined in a bigger domain. Using the relation 

Trln{I - uM) = lndet{I - uM) (4.53) 

which is vahd for det{l — uM) [18], we get 

^TTT-T^ / / dedr)lndet(I - uM) (4.54) 
2 (27r)"' Jo Jo 

The determinant can be easily computed and one finds that 

det(l - uM) = (m^ + if - 2u{l - u^){cose + cost;) (4.55) 

Taking the logarithm on both sides of relation (3.24) gives 

^^^1^ = ln2 + 2(1 - hln{coshK) + In + Wp{u)] (4.56) 
i V iv p 

or, using (4.45-47), (4.52,4.54), 



lnZ{u) 



ln2 + 2(1 - —)ln(coshK) + At /^^ /^^ dedr]lndet(l - uM) (4.57) 
N 87r"' Jo Jo 
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Using (4.55) with u — tanhK and the relations 

1 



gives 



u = -sinh{2K){l - u^) (4.58) 
[l + u^ = cosh''{2K){l-u^f (4.59) 



det{l - uM) = cosh~^2K {cosh2Ky - {sinh2K){cosri + cose) (4.61) 



Onsager's formula follows from (4.57) after taking the limit N — >■ oo: 

--^^ln2 + ^j I dedr]ln\(cosh2Kf - sinh2K( COST] + cose)] (4.62) 
kbT 2tt^ Jo Jo ^ -i 

where / is the free energy per site in the thermodynamical limit ( see (2.5)). 

The integral in (4.62) can not be evaluated in terms of simple functions. The 
derivatives of the integral however can be expressed in terms of elliptic functions 
[2,15,19]. 

Set 2k = tanh2K / cosh2K . Then, 



/ 



ln2 + ln(cosh2K) + -— r / / ded7]ln [1 - 2k(cos7] + cose)] (4.63) 

277"^ Jo Jo 



Expanding the logarithm in powers of k it follows that 



^ + in(co./^2X) - £ (1^)' fc'" (4.64) 

The series converges for |2A;(cose + cos77)| < 4|/c| < 1. For A; > (J > 0) and at A; = 
1/4, that is, at the critical value K — Kc { or temperature Tc — 2Jk'^^lrr^{-\/2 -\- 1) 
) given by 

2senh2K^ = cosh^2K^ (4.65) 

it diverges. ( Similarly, for J < which implies A; < and divergence at Tc = 

2JksHn-\V2 - 1) ). 

The internal energy U is given by 



.2d f 



U = -kT' — ^ (4.66) 
Prom (4.63), 

U = -Jcoth{2K) [l + {sinh^2K - l)/(2i^)] (4.67) 
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where 

, , 1 /■'^ dedr] 

^ ■'^n^Jo cosh^{2K) - sinh{2K){cose + cos?]) ^ ' ' 

By performing one of the integrals its found that 



U = -Jcoth{2K) 



1 + {2tanh'^2K -1)-F{ki) 

TT 



(4.69) 



where ki — Ak and F{ki) is the complete elliptic integral of the first kind defined by 

F(ki)= r^\l-klsen^9)-U9 (4.70) 
Jo 

The elliptic function F ( see Ref. [19] ) has the property that 

F ^ ln[4{l - kf)]-^/^ (4.71) 

as /ci — > 1~. So it diverges logarithmically at /ci = 1, or at the value Kc given by 
(4.65). 

In relation (4.69) for U the function F is multiphed by {2tanh'^2K — 1) which is 
zero at the critical point K^. Indeed, from the identity cosh^x = 1 + sinh'^x relation 
(4.65) imphes that sinh{2Kc) = 1. Using this and (4.65), tanh^2Kc — 1/2 follows. 
So the function U is continuous at Kc- 

The specific heat can be computed from the definition 

(4.72) 

It is given by 
OK 

C = —{Kcoth2K)'^{2F{ki) - 2E{ki) + 2{tanh?2K - l)G{ki)} (4.73) 

TT 



where 



G{ki) = ^ + {2tanh^2K - l)F{k{) (4.74) 
and E{ki) is the complete elliptic integral of the second kind, defined by 

E{ki) = r^\l - klsen'^e)Ue (4.75) 

which is well defined at ki — 1. Prom the exact result (4.73) it follows that the 
specific heat is logarithmically divergent at the critical point. 
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